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(57) ABSTRACT 

Method and system for estimating a growth stage of an alga in 
an ambient fluid. Measured light beam absorption or reflec- 
tion values through or from the alga and through an ambient 
fluid, in each of two or more wavelength sub-ranges, are 
compared with reference light beam absorption values for 
corresponding wavelength sub -ranges for in each alga growth 
stage to determine (1 ) which alga growth stage, if any, is more 
likely and (2) whether estimated lipid content of the alga is 
increasing or has peaked. Alga growth is preferably termi- 
nated when lipid content has approximately reached a maxi- 
mum value. 


16 Claims, 11 Drawing Sheets 
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When (1) e (g = gO) for a particular growth stage, g = gO is less than or equal to 
e(g) for any other growth stage g and (2) e (gO) < e(thr), interpret these 
conditions as indicating that the alga is more likely in the reference growth stage 
g = gO than in any other growth stage 
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ESTIMATION OF ALGA GROWTH STAGE 
AND LIPID CONTENT GROWTH RATE 

ORIGIN OF THE INVENTION 

5 

This invention described herein was made in the perfor- 
mance of work under a NASA contract and by an employee of 
the United States Government and is subject to the provisions 
of Public Law 96-517 (35 U.S.C. §202) and maybe manufac- 
tured and used by or for the Government for governmental 
purposes without the payment of any royalties thereon or 
therefore. 

FIELD OF THE INVENTION 

15 

This invention relates to methods of estimation of alga 
growth stage and lipid content growth rate for alga in marine 
water, in fresh water and in any other fluid, using light absorp- 
tion or light reflection (collectively referred to as “light modi- 2 o 
fication”) in selected wavelength ranges. 

BACKGROUND OF THE INVENTION 

In many approaches to production of fuels, nutraceuticals 25 
and other useful products by growth and conversion of algae 
products, a growth stage of alga in a controlled environment 
must be evaluated, for purposes of estimation of lipid contant 
and for other metrics. It is often very inaccurate to estimate 
alga growth stage based only on time of growth, because of 30 
variations in important environmental parameters during dif- 
ferent growth cycles. 

What is needed is a more accurate approach that does not 
depend solely on time of growth and that implicitly factors in 
the variable environmental parameters and provides a useful 35 
error estimate for alga growth stage. Preferably, the approach 
should implicitly or explicitly take account of the environ- 
mental history of growth of the alga. 

SUMMARY OF THE INVENTION 40 

These needs are met by the invention, which provides a 
method using selective modification of light in different 
wavelength ranges to estimate (i) alga growth stage and (ii) 
alga growth rate in an ambient fluid (e.g., air, vacuum, fresh 45 
water, marine water and/or brine). In a first embodiment, 
absorption of light is measured for a beam, having a specified 
light intensity in each of two or more specified narrow wave- 
length ranges of light, Xl n ^X^X2 n . (n=l, 2, . . . , N; Ni^3). 
Optionally, light absorption is compensated for absorption of 50 
light in the same wavelength range by the ambient medium 
(vacuum or a specified fluid). The net absorption of light in N 
wavelength ranges (Ni^3) is compared with a reference set of 
net absorption values for the alga for different growth stages. 

An algorithm is applied to differences between measured net 55 
absorption values and reference net absorption values to esti- 
mate a growth stage (one of a discrete sequence of growth 
stages). Compensation for light reflection from a liquid (ab- 
sent alga) is similar. 

In a second embodiment, lipid content of the alga is mea- 60 
sured or estimated at each of a selected set of growth stages. 

In some algae, the lipid content increases to a maximum, and 
the lipid content thereafter reaches a plateau. In some other 
algae, the lipid content increases to a maximum and thereafter 
decreases to a lower value, indicating that the alga growth 65 
process should be terminated at some time in order to accu- 
mulate the largest lipid content from that alga. 
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In a third embodiment, the estimated growth stage is cor- 
related with a time variable, t-tO, where tO is an estimated 
time for initiation of growth the alga under specified condi- 
tions. A temporal rate of alga growth from one stage to the 
next stage is estimated, and this rate is correlated with accu- 
mulated time, t-tO, to allow an estimate of sensitivity of 
accumulated lipid content with a time for cut-off of alga 
growth. 

In a fourth embodiment, one or more relevant environmen- 
tal parameters (light intensity, light wavelength, temperature, 
iron content, etc.) is varied for the alga, and a time At(hours) 
required for the algae to progress from an initial stage to a 
specified “end stage” is determined, based on variation of the 
parameters one at a time. The inverse l/At(hours) is taken to 
represent an average growth rate of the algae. A particular 
combination of environmental parameter values is identified 
for which the average growth rate is greatest (l/At(hours) is 
smallest). 

BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. 1 schematically illustrates an experimental set-up for 
practicing an embodiment of the invention. 

FIG. 2 is a flow chart of a procedure for practicing the 
invention. 

FIGS. 3-8 are graphical illustrations of alga absorption 
versus wavelength for selected algae. 

FIGS. 9-10 are graphical illustrations comparing alga 
specitic mass and alga lipid content versus alga growth time. 

DESCRIPTION OF BEST MODES OF THE 
INVENTION 

FIG. 1 schematically illustrates a system 11 for practicing 
the invention. A specified alga AL is suspended in a chamber 
12 in an ambient fluid medium 13 (vacuum, air, fresh water, 
marine water, a specified gas or vapor at a specified pressure, 
etc.), where the chamber 12 has first and second windows (or, 
optionally, a single wraparound window), 14-1 and 14-2, that 
face each other and that are substantially transparent through- 
out a range of wavelengths, X(lower)^X^X(upper), that 
includes the wavelengths of interest. A light beam 15, pref- 
erably filtered to transmit light in only one of a sequence of N 
wavelength ranges, X\ n ^k^k2 n . (n=l, . . . , N; N^3), illu- 
minates the alga AL through the first widow 14-1. Light not 
modified by the alga and not absorbed (or reflected) by the 
ambient medium within the chamber 12 passes through the 
second window 14-2, or a second portion of the first window, 
and modified light is received by a light measurement mecha- 
nism 16, which sequentially or simultaneously measures light 
intensity I(n) received in the wavelength ranges X l n ^X^X2 n . 
(n=l, . . . , N; Ni=3) 

A computer 17 associated with the light measurement 
mechanism 16 receives the measured (modified) light inten- 
sity value I(n) for each of N wavelength ranges, Xl n ^X^X2 n , 
and estimates light modification in each wavelength range of 
interest. Optionally, the computer 17 applies a procedure to 
compensate for light intensity modified by the ambient 
medium (with alga absent) and by the window material in 
each wavelength range of interest, Xl n ^X^X2 n . 

For example, if a simple Beers’ law exponential absorption 
model is applied, the total light absorption is a sum of the 
exponential absorption factors, a(alga;n) and a(ambient;n) 
and a(window;n), for the alga and for the ambient medium 
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and window(s), respectively; and the net alga transmission is 
then expressible as 

/(net;«)=/ 0 (K)exp{a(ambient;K)+a(wmdow;K)}exp{- 

a(alga;«)}, (1) 

with a(alga;n) replaced by a(alga refl;n) for light reflection in 
Eq. (1). The factor exp {a(ambient;n)+a(window;n)} is esti- 
mated in a separate experiment and is made available as a 
reference value for different algae. The computer 17 com- 
putes the net transmission I(net;n)* and estimates net alga 
modification in the wavelength range, W n ^k^'k2 n . Option- 
ally, the computer 17 can estimate a volume or opaqueness 
factor for the alga (displacement of fluid by the alga) within 
the chamber and can further modify the alga modification 
factor a(alga;n) to account for absence of the ambient fluid 
wherever the alga is present. 

A suitable wavelength range width for measurement of 
alga absorption may be X2 w -Xl w =5-50 nm, depending upon 
the average wavelength, {X2„+Xl„}//2. Where the absorption 
factor a(alga;n) is changing rapidly with the index n, the 
wavelength range width should be smaller; use of non-uni - 
form widths is preferred here. 

A sequence of net absorption values {a(alga;n)}„ for each 
of the wavelength ranges of interest is compared with a ref- 
erence set of values {a(ref;n;alga growth stage g)} K to esti- 
mate an error, between reference value and measured value, 
for each of G stages of growth, numbered g=l, . . . , G (G^2). 
FIGS. 3-8 illustrate a representative group of N net light 
modification values (N=3) for different growth stages of an 
alga, Chlorella vulgaris, Chlorella (mixed species) and 
Scenedesmus. 

One estimate of the associated error is 


* (2) 
s(g) = \ w n \a{net\ n ) - a{ref\ n\ growth stage = 1, ... , G ), 

n = 1 

where {w w } w is a sequence of selected non-negative weight 
values and p is a selected positive number, for example, p=l, 
p=2 or p=selected rational or irrational number. The growth 
stage error values e(g) are compared with each other and with 
a threshold error value e(thr). Where a particular growth stage 
error value, €(g=gO), satisfies the conditions 


e(^)=mm{e(g=l), . . . , e(g=G)}, 

(3A) 

e(g0)^e(thr), 

(3B) 


this is interpreted as indicating that the alga growth stage 
g=gO is the most likely growth stage, based on the N mea- 
surements of alga net absorption. 

FIG. 2 is a flow chart of a procedure for practicing an 
embodiment of the invention. In step 21, an alga, having a 
growth stage that may be unknown, is placed in a chamber 
that also contains an ambient fluid (liquid or gas or vacuum), 
the chamber having one or more windows that are substan- 
tially transparent to light in a selected wavelength range R. 

In step 22. N light beams, having initial light beam inten- 
sities I 0 (n) (n=l, . . . , N; N^2) within the wavelength sub- 
ranges, respectively, within the range R, are passed through 
the one or more chamber windows, through the ambient fluid 
and through the alga, and at least one of the ambient fluid and 
the alga is allowed to absorb a portion of at least one of the N 
light beams, to produce N modified light beam intensities I(n) 
that issue from the chamber. 

In step 23, the N modified light beam intensities are 
received and measured or estimated. 
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In step 24, the N modified light beam intensities I(n) are 
compared with the corresponding N initial light beam inten- 
sities I 0 (n) 

In step 25, N reference light beam intensities I(ref;n;g) are 
5 provided for the N wavelength ranges for each reference 
growth stage g (g=l, . . . , G;G^2). 

In step 26, an error value €(g) is computed, based upon 
differences between the N modified light beam intensities and 
the corresponding N reference light beam intensities for each 
to of the reference alga growth stages, g=l, . . . , G, 

When (1) the error value €(g=g0) for a particular growth 
stage, g=g0, is less than or equal to the error value €(g) for any 
other value of g and (2) the error value e(g0) is no greater than 
a selected threshold error value €(thr), these conditions are 
15 interpreted, in step 27, as indicating that the alga is more 
likely in the reference growth stage g=g0 than in any other 
reference growth stage. 

When (3) the error value e(g0) is greater than the selected 
threshold error value €(thr), these conditions are interpreted, 
20 in step 28, as indicating that it cannot be determined, from 
these conditions alone, whether the alga is more likely in a 
particular one of the reference alga growth stages. 

In step 29 (optional), the lipid content LC(t m ) of the alga is 
estimated for each of a sequence of absorption measurement 
25 times, t=t w (m=l, . . . , M 0 , and the estimated value LC(t m ) is 
associated with the estimated alga growth stage, g=g0, deter- 
mined in the preceding step. In step 30 (optional), the system 
determines whether the estimated lipid content LC(t m ) 
increases (strictly), expressed as LC(t m+1 )i=LC(t m ). Where 
30 the answer to the query in step 30 is “yes ,” the system inter- 
prets this response as indicating that the estimated lipid con- 
tent does not (yet) decrease with increasing time, t=t m , and 
alga growth is allowed to continue, in step 31 (optional). 

Where the answer to the query in step 30 is “no” so that the 
35 estimated lipid content is no longer increasing (strictly) 
monotonically, the system interprets this response as indicat- 
ing that estimated lipid content may decrease with time 
beyond a certain growth stage. One consequence of this last 
interpretation is that alga growth should be terminated, in step 
40 3 2 (optional), before growth continues, in order to maximize 
lipid content that may be extracted from the alga. 

Certain algae, when fully grown, have a relatively high 
lipid content and can be used as a feedstock for extraction of 
biofuels, diesel fuels, nutraceuticals, pharmaceuticals, fertil- 
45 izer, animal feed and other useful products. Algae to which 
the invention has been applied to evaluate alga growth stage 
and alga lipid content include Chlorella vulgaris, Chlorella 
(mix) and Scenedusmus. 

FIGS. 3 and 4 are graphical illustrations of measured 
50 absorption a(alga;n) in Chlorella vulgaris versus wavelength 
X(nm) for a wavelength range of about 300 nm^X^750 nm at 
four days growth (FIG. 3) and at nine days growth (FIG. 4). 
The following differences are observed between FIGS. 3 and 
4: (1) a local maximum in absorption value at X«370 nm for 
55 four days growth is larger than the corresponding peak for 
nine days growth; (2) a local minimum in absorption at X«390 
nm for four days growth has a larger value than the corre- 
sponding absorption minimum value for nine days growth; 
(3) the maximum absorption value (at X«400 nm) for nine 
60 days growth is larger than the corresponding absorption peak 
for four days growth; (4) a local minimum absorption value 
(at X«420 nm) has a higher value for nine days growth than 
the corresponding local minimum absorption value for four 
days growth; (5) a local maximum absorption value (at X«700 
65 nm) for four days growth has a higher value than the corre- 
sponding local maximum absorption value for nine days 
growth. 
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FIGS. 5 and 6 are graphical illustrations of measured 
absorption a(alga;n) in Chlorella (mixed species) versus 
wavelength X(nm) for a wavelength range of about 300 
nm^X^750 nni at growth stage four days growth (FIGS. 5) 
and 16 days growth (FIG. 6). The following differences are 
observed between FIGS. 5 and 6: (1) a local absorption maxi- 
mum (at X«380 nm) for 16 days growth is more pronounced 
than the corresponding peak at four days growth; (2) a sub- 
sequent minimum absorption value (at X=395 nm) for 1 6 days 
growth is larger than a corresponding minimum absorption 
value for four days growth; (3) the wavelength for the mini- 
mum value in (2) for 16 days growth is shifted relative to the 
wavelength for the corresponding minimum value for four 
days growth; (4) the maximum absorption peak (at X«470 
nm) for 16 days growth is larger than the corresponding 
maximum absorption peak for four days growth; (5) a subse- 
quent local absorption maximum value (at X«520 nm) for 1 6 
gays growth is larger than the corresponding peak for four 
days growth; (6) a local minimum value for absorption (at 
X«680 nm) for four days growth has disappeared at 16 days 
growth; and (7) a local maximum absorption peak (at X«710 
nm) for four days growth is larger than the corresponding 
local maximum absorption peak for 16 days growth. 

FIGS. 7 and 8 are graphical illustrations of measured 
absorption a(alga;n) in Scenedesmus versus wavelength 
X(nm) for a wavelength range of about 300 nm ^X.^750 nm at 
growth stage seven days growth (FIG. 7) and 1 5 days growth 
(FIG. 8). The following differences are observed between 
FIGS. 7 and 8: (1) a local maximum absorption value (at 
X«380 nm) for 1 5 days growth is larger than the correspond- 
ing local maximum value for seven days growth; (2) a local 
minimum absorption value (at X«395 nm) for 1 5 days growth 
is larger than the corresponding local minimum value for 
seven days growth; (3) the maximum absorption value (at 
X«480 nm) for 1 5 days growth is larger than the correspond- 
ing maximum absorption value for seven days growth; (4) a 
local minimum absorption value (at X«500 nm) for 15 days 
growth is larger than the corresponding minimum value for 
seven days growth; (5) a local maximum absorption value (at 
X«500 nm) for 1 5 days growth is larger than the correspond- 
ing maximum value for seven days growth; and (6) a local 
maximum absorption value (at X«690 nm) for seven days 
growth is larger than the corresponding maximum value for 
1 5 days growth, where this local maximum value has nearly 
disappeared. 

One or more of five factors contributes to absorption graph 
differences between two or more distinct alga growth stages: 
(i) differences in (local or absolute) absorption maximum 
values; (ii) differences in wavelength value for the (local or 
absolute) absorption maximum values; (iii) differences in 
local absorption minimum values; (iv) differences in wave- 
length value for the local absorption minimum values; and (v) 
appearance of a prominent feature (e.g., maximum or mini- 
mum absorption value) at one growth stage that is absent at 
another growth stage. Reflection graph differences between 
two alga growth stages arise from: (i) appearance or disap- 
pearance of a prominent reflection feature at a growth stage 
and (ii) reflectivity differences at different growth stages. 

FIGS. 9 and 10 graphically illustrate alga mass (gms/liter) 
and percent lipid content for Chlorella vulgaris , for a bag 
filled with City of Sunnyvale primary effluent plus C0 2 (FIG. 
9) and for a bag filled with BG-1 1 media plus C0 2 (FIG. 10), 
versus number of days of alga growth (1-14). With reference 
to FIG. 9, the alga specific mass increases to a value of about 
1 5 gm/liter at about seven days growth and thereafter plateaus 
or decreases slightly, while lipid content of the alga continues 
to increase to about 22 percent at about ten days growth, and 


6 

thereafter decreases sharply and plateaus at a lower value of 
about 1 2 percent on day 1 4. These responses indicate that, for 
Chlorella vulgaris , in the Sunnyvale effluent, the alga should 
be harvested shortly after the alga first achieves maximum 
5 growth, in order to obtain maximum lipid content. 

With reference to FIG. 10, the alga specific mass increases 
through day 12; the lipid percent increases until day 10 and 
decreases sharply beyond day 1 0, when the alga specific mass 
is still increasing. These responses indicate that, for Chlorella 
10 vulgaris , in the BG-11 media, the alga should be harvested 
before or shortly after the alga first achieves maximum 
growth, in order to obtain maximum lipid content. Delaying 
the harvest beyond about nine or ten days growth will result in 
reduced total lipid content. The differences inn the lipid con- 
1 5 tent response in FIGS . 9 and 1 0 indicate that lipid growth may 
depend upon the particular alga and upon the alga growth 
medium, 

From FIGS. 9 and 10, one notes that the rate of alga growth 
varies substantially with its present “age,” and accurate pre- 
20 diction of ultimate (saturated) specific mass of the alga is 
probably not possible. However, one can estimate a present 
rate of alga growth, using two or three rates at preceding times 
as a predictor. When the estimated present rate of growth is 
below a specified threshold (e.g., 2-4 gm/liter/day), one can 
25 conclude that alga growth saturation is near and treat the 
current state of the alga as an end state. 

Estimated lipid content associated with the alga growth at 
an alga harvest time can also be optimized by the following 
procedure. The rate of alga growth at each of the sequence of 
30 measurement times, t=t m , is estimated as follows. For three 
measurement times, t=t m _ l5 t m , and t m+1 , with t m _ 1 <t 1 <t m+1 
the quadratic function 

Q(t; t m - 1, t m , r„ + i,)=iC(f M _ 1 ) (f-UH + i)/(^-i-U 

35 (^-^ + i)+^C(r ra+1 )(r-r m _ 1 )(r-0/(^_i-^_i) 

( 4 ) 

reproduces the lipid content values LC(t m+x ) (x=-l, 0, +1) at 
the respective time values t=t m _ x . The time derivative of Q(t; 

t m _ i, t m , t m+lJ ) at the intermediate time value, t=t m ), is 
40 

(a Q/dt^^LCit^) (^-^ +1 )/(^_ 1 -^)F ra _ 1 - 
t m +i)+LC(t m ) 

t m+i)+LC(t m+l ) ( 5 ) 

When the lipid content growth rate {3Q/3t}^ m is negative, is 
45 zero, or is positive but small relative to preceding values of the 
lipid content growth rate, alga growth should be terminated 
and the alga (and associated lipid content) should be har- 
vested, at a time contemporaneous with the alga termination 
time t(term). A computer can be programmed to compute the 
50 quantity {3Q/3t}^ m at each of the measurement times, t=t m , 
and to determine when to terminate alga growth. 

What is claimed is: 

1. A method for estimating a stage of growth of a selected 
algae, the method comprising: 

55 placing an algae, having a growth stage that may be 
unknown, in a chamber that also contains an ambient 
fluid, the chamber having one or more windows that are 
at least partly transparent to light in a selected wave- 
length range R; 

60 passing N light beams, numbered n=l, . . . , N (N^2), 
having initial light beam intensities I 0 (n) and having 
wavelength sub-ranges 'k\ n ^k^K2 n ., respectively, 
within the range R, through the one or more chamber 
windows, through the ambient fluid and through the 
65 algae, and allowing at least one of the ambient fluid, the 

one or more chamber windows and the algae to absorb a 
portion of at least one of the N light beams, to produce 
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modified light beams, numbered n=l, . . . , N, in the 
respective wavelength sub-ranges, that have passed 
through the chamber; 

receiving the modified N light beams and estimating or 
measuring modified light beam intensities, I(n), num- 
bered n=l, . . . , N, for the respective N modified light 
beams; 

providing N reference light beam intensities I(ref;n;g), 
numbered n=l , . . . , N, for each reference growth stage, 
numbered g=l, . . . , G (Gi=2) for the selected algae; 

comparing the N modified light beam intensities I(n) with 
the corresponding N reference initial light beam inten- 
sities I(ref;n;g) for each reference growth stage g of the 
algae; 

computing an error value €(g), based upon differences 
between the N modified light beam intensities I(n) and 
the respective N reference light beam intensities I(ref;n; 
g), for each growth stage g of the algae; and 

when, for a selected growth stage number gO, (1) the error 
value €(g) for at least one growth stage, g-gO, is no 
greater than a value €(g'), for any other growth stage g’, 
and (2) €(g0) is no greater than a selected error threshold 
value €(thr), interpreting these conditions as indicating 
that the algae is most likely in the growth stage gO. 

2 . The method of claim 1 , further comprising: 

when, for said selected algae growth stage number gO, (3) 
€(g) is greater than €(thr) for said selected error threshold 
value, interpreting this condition as indicating that it 
cannot be determined, from this condition alone, 
whether said algae is most likely in said growth stage gO. 

3 . The method of claim 1 , further comprising choosing at 
least one of said error value e(g) to be expressed as 


N 

£{g) = z w n \a{n\ net ) - a{ref ; n\ growth stage < g r )| p , 


where g refers to said growth stage of said algae, a(n;net) is a 
light beam absorption factor corresponding to said modified 
light beam intensity I(n), a(ref;n; growth stage g) is a refer- 
ence light beam absorption factor corresponding to said ref- 
erence light beam intensity I(ref;n;g) for said algae growth 
stage g, w w is a selected non-negative weight value, and p is a 
selected positive number. 

4 . The method of claim 1 , further comprising choosing said 
ambient fluid to include at least one of air, a vacuum, fresh 
water and marine water. 

5 . The method of claim 1 , further comprising: 
providing an estimate LC(g’) of lipid content in said algae 

at each of a sequence of growth stages g—1, . . . , G’; 
estimating a growth stage, g ,= g'(max) at which the lipid 
content LC(g') is maximized; and 
when said most likely growth stage gO has a value at least 
equal to g'(max), terminating a growth process for said 
algae. 

6. The method of claim 1, further comprising: 
providing an estimate LC of lipid content in said algae at 

each of a sequence of growth stages g', corresponding to 
times t=t g ., (g'=l, . . . , G'; G'g3}; 
estimating a time rate of change, 3LC/3t, of the lipid con- 
tent LC for at least one selected time, t=t^,, t=t.^, +1 and 
t=V +2 ; and 

when the estimated time rate of change of the lipid content 
for the selected time, t=t^,, is either negative or 0 or has 
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a small positive value, terminating growth of said algae 
at a time contemporaneous with the time t=t 

7. The method of claim 6, wherein said process of estimat- 
ing said time rate of change 3LC/3t comprises: 

5 providing a sequence {LC(t m )} m of measurements or esti- 
mates of lipid content of a selected alga at a sequence of 
times, t=t m , and 

for at least three measurement times, t=t m _ 1? t m , t m+1 , with 
estimating said lipid content growth rate 
10 at a time, t=t m , as 

{dQ/dtj^LC^) 

^ + i) + ZC(r m )(2r ra -r ra _i-r ra+1 )/(^-^_ 1 )(^- 

^ + i)+ZCti ra+1 )(r-r m _ 1 )/(r w _-^_ 1 )(r ra+1 -0- 

8. The method of claim 1, wherein said process of provid- 
ing N reference light beam intensities I(ref;n;g) for at least 
one of said algae growth stages g for said selected algae 
comprises 

20 providing an exponential absorption value a(ambient;n) 
for a thickness of said ambient fluid in which said algae 
is immersed; 

providing an exponential absorption value a(window;n) 
for said chamber window or chamber windows through 
25 which at least one of said light beams passes; 

estimating an overall light beam attenuation as exp {-a 
(algae;n;g)}, and 

estimating said reference light beam intensity by a relation 
I(net;n;g)=I 0 (n) exp{a(ambient ;n)+a(window ;n) } 

30 exp{-a(algae;n)}. 

9. A method for estimating a stage of growth of a selected 
algae, the method comprising: 

placing an algae, having a growth stage that may be 
unknown, in a chamber that also contains an ambient 
35 fluid, the chamber having one or more windows that are 
at least partly transparent to light in a selected wave- 
length range R; 

passing N light beams, numbered n=l, . . . , N (Ni=2), 
having initial light beam intensities I 0 (n) and having 
40 wavelength sub-ranges X1 W ^X^X2 W , respectively, 

within the range R, through the one or more chamber 
windows, through the ambient fluid, allowing the algae 
to reflect the light as a reflected light beam, and allowing 
at least one of the ambient fluid, the one or more chamber 
45 windows and the algae to reflect a portion of at least one 

of the N light beams, to produce modified light beams, 
numbered n=l, . . . , N, in the respective wavelength 
sub-ranges, that have been reflected; 

receiving the modified N light beams and estimating or 
50 measuring modified light beam intensities, I(n), num- 

bered n=l, . . . , N, for the respective N modified light 
beams; 

providing N reference light beam intensities I(ref;n;g), 
numbered n=l, . . . , N, for each reference growth stage, 
55 numbered g=l, . . . , G (G^2) for the selected algae; 

comparing the N modified light beam intensities I(n) with 
the corresponding N reference initial light beam inten- 
sities I(ref;n;g) for each reference growth stage g of the 
alga; 

60 computing an error value €(g), based upon differences 
between the N modified light beam intensities I(n) and 
the respective N reference light beam intensities I(ref;n; 
g), for each growth stage g of the algae; and 

when, for a selected growth stage number gO, (1) the error 
65 value €(g) for at least one growth stage, g-gO, is no 
greater than a value e(g'), for any other growth stage g T , 
and (2) €(g0) is no greater than a selected error threshold 
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value €(thr), interpreting these conditions as indicating 
that the algae is most likely in the growth stage gO. 

10. The method of claim 9, fiirther comprising: 

when, for said selected growth stage number gO, (3) €(g) is 
greater than €(thr) for said selected error threshold value, 5 
interpreting this condition as indicating that it cannot be 
determined, from this condition alone, whether said 
algae is most likely in said growth stage gO. 

11. The method of claim 9, further comprising choosing at 
least one of said error value €(g) to be expressed as 


N 

s(g ) = ^ w n \a{n\ net ) - a{ref\ n\ growth stage g)| p , 

n= l 

15 

where g refers to said growth stage of said algae, a(n;net) is a 
light beam modification factor corresponding to said modi- 
fied light beam intensity I(n), a(ref;n; growth stage g) is a 
reference light beam absorption factor corresponding to said 
reference light beam intensity I(ref;n;g) for said algae growth 20 
stage g, w„ is a selected non-negative weight value, and p is a 
selected positive number. 

12. The method of claim 9, further comprising choosing 

said ambient fluid to include at least one of air, a vacuum, 
fresh water and marine water. 25 

13. The method of claim 9, fiirther comprising: 
providing an estimate LC(g') of lipid content in said algae 

at each of a sequence of growth stages g— 1, . . . , G’; 
estimating a growth stage, g ,= g'(max) at which the lipid 
content LC(g') is approximately maximized; and 30 

when said most likely growth stage gO has a value at least 
equal to g'(max), terminating a growth process for said 
algae. 

14. The method of claim 9, further comprising: 
providing an estimate LC of lipid content in said algae at 35 

each of a sequence of growth stages g T , corresponding to 
times t t A ,,, (g'=l, . . . , G'; G'g3}; 
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estimating a time rate of change, 3LC/3t, of the lipid con- 
tent LC for at least one selected time, t=t * t=t, , +1 and 
t=V +2 ; and 

when the estimated time rate of change of the lipid content 
for the selected time, t=t^,, is either negative or 0 or has 
a small positive value, terminating growth of said algae 
at a time contemporaneous with the time t=t^,. 

15. The method of claim 14, wherein said process of esti- 
mating said time rate of change aLuat comprises: 

providing a sequence {LC(t m )} m of measurements or esti- 
mates of lipid content of a selected alga at a sequence of 
times, t=t m , and 

for at least three measurement times, t=t m _ 1; t m , t m+1 , with 
t m _ 1 <t m <t m+1 , estimating said lipid content growth rate 
at a time, t=t m , as 

^ + i)+2:c(^ +1 )(r-r w _ 1 )/(? ra _-^_i)(^ + i-0- (5). 

16. The method of claim 9, wherein said process of pro- 
viding N reference light beam intensities I(ref;n;g) for at least 
one of said algae growth stages g for said selected algae 
comprises 

providing an exponential absorption value a(ambient;n) 
for a thickness of said ambient fluid in which said algae 
is immersed; 

providing an exponential absorption value a(window;n) 
for said chamber window or chamber windows through 
which at least one of said light beams passes; 

estimating an overall light beam attenuation as exp {-a 
(algae;n;g)}, and 

estimating said reference light beam intensity by a relation 
I(net;n;g)=I 0 (n) exp { a(amb ient ;n)+a(window ;n) } 

exp{-a(algae;n)}. 





